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Abstract 

We  summarize  progress  in  our  analysis  of  proposed  communication  schemes  for  hypersonic  ve¬ 
hicles  in  flight.  Investigations  into  the  electron-acoustic  wave  (EAW)  communication  scheme  have 
yielded  a  dispersion  analysis  indicating  that  EAW  modes  can  be  generated  and  propagate  in  the 
plasma  layer  [D.  V.  Rose  et  al,  Voss  Sci.  Rep.  No.  VSL-0737,  Nov.  2007].  These  modes  can  be 
coupled  to  electromagnetic  waves  at  the  plasma  layer  boundary  for  the  assumption  of  a  sharp  den¬ 
sity  boundary.  A  detailed  analysis  of  the  wave  transformation  coefficients  at  this  boundary  layer 
has  been  completed  [D.  V.  Rose  et  al,  Voss  Sci.  Rep.  No.  VSL-0832,  Dec.  2008].  Here  we  report  on 
new  particle-in-cell  simulations  of  EAW  generation  and  propagation  in  a  two-temperature  plasma. 
The  simulation  model  is  being  developed  to  test  the  stability  and  the  transformation  efficiencies  of 
the  waves  in  plasma  conditions  anticipated  for  hypersonic  vehicles  in  sustained  upper  atmospheric 
flight.  In  addition,  we  summarize  progress  in  our  analysis  of  the  ReComm  scheme  [M.  Keidar  et  al., 
J.  Spacecraft  and  Rockets  45,  445  (2008)],  focusing  on  extending  our  previous  ID  computational 
model  [C.  Thoma  and  D.  V.  Rose,  Voss  Scientific  Report  No.  VSL-0816,  July  2008]  to  longer  time 
and  length  scales. 
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I.  INTRODUCTION 


Hypersonic  vehicles  traveling  in  the  upper  atmosphere  at  speeds  greater  than  Mach  10 
(7,000  mph  at  33  km  altitude)  generate  plasma  that  disrupt  or  prevent  communications  over 
conventional  radio-frequency  channels.  This  phenomenon  manifests  itself  most  famously  as 
the  well-known  communications  “blackout”  period  during  space  vehicle  reentry  into  the 
atmosphere  [1-3],  and  has  been  studied  (off  and  on)  since  the  Gemini  and  Apollo  space 
programs.  For  sustained  hypersonic  flight  in  the  atmosphere,  the  communications  blackout 
will  persist  for  almost  the  entire  flight. 

The  plasma  boundary  layer  is  formed  around  a  hypersonic  vehicle  in  the  bow  shock  and 
to  some  extent  by  the  collision  of  neutral  gas  particles  with  material  ablated  from  the  surface 
of  the  hypersonic  vehicle  itself.  In  the  shock,  the  air  is  heated  to  such  a  high  temperature 
that  ionizing  collisions  between  neutral  particles  occur.  Estimates  of  electron  density,  Up 
(~  10^^  cm“^),  temperature  (~  0.5  eV)  and  electron-neutral  collision  frequencies  (2.5  x  10®- 
42.5  X  10^^  s“^)  have  been  obtained  for  high-Mach-number  reentry  vehicles  at  altitudes 
between  100,000-150,000  ft  (30-45  km)  [4,  5].  The  thickness  of  the  plasma  layer  (defined  as 
the  region  where  the  electron  density  exceeds  ~  10®  cm“^)  associated  with  these  parameters 
is  on  the  order  of  6  cm,  with  90%  of  the  peak  density  out  to  a  width  of  about  2  cm.  The  free 
electrons  in  the  plasma  layer  attenuate  radio- frequency  (RF)  waves  both  through  reflection 
and  resistive  absorption,  and  generate  RF  noise.  At  sufficiently  large  antenna  power  levels, 
the  RF  itself  can  cause  further  ionization  of  the  air  to  occur. 

Acceptably  low  attenuation  of  electromagnetic  (EM)  wave  propagation  in  a  collisionless 
plasma  is  limited  to  frequencies  above  the  cutoff  value  fc  =  Upe/2n,  where  Upe  is  the  electron 
plasma  frequency  (see,  for  example.  Ref.  [6]).  For  Up  =  10^^  cm“^,  the  fc  ~  10  GHz,  which 
is  at  the  high  end  of  the  radio  wave  frequency  spectrum  and  well  into  the  conventional 
microwave  portion  of  the  EM  spectrum.  In  addition,  this  plasma  cutoff  frequency  is  above 
the  spectrum  allocated  (in  some  cases  by  international  treaty)  for  flight  test  and  evaluation 
telemetry  and  above  frequencies  supported  by  existing  infrastructure.  We  note  that  Global 
Positioning  System  signals  are  broadcast  in  the  L-band  region  of  the  spectrum,  and  con¬ 
sequently  are  also  heavily  attenuated  by  the  plasma  sheath.  As  a  result,  communication 
approaches  within  the  standard  telemetry  spectrum  are  inadequate  for  hypersonic  vehicles 
using  conventional  techniques  (see,  for  example.  Refs.  [6-9]). 
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A  theoretical  research  program  is  underway  to  assess  concepts  for  real-time  telemetric 
communications  through  plasma  boundary  layers  surrounding  hypersonic  vehicles  in  flight. 
This  research  includes  detailed  analytic  analysis  and  computational  modeling  of  several 
proposed  communication  schemes.  In  our  previous  reports  [10-14],  we  summarized  our 
research  progress  on  fonr  such  schemes:  a  nonlinear  three-wave  interaction  techniqne  [15- 
17],  an  electron-acoustic  wave  (EAW)  propagation  scheme,  a  magnetic  window  scheme  for 
EM  wave  propagation  throngh  a  collisional  plasma  slab  [12],  and  a  whistler  wave  excitation 
mode  [11].  Here  we  snmmarize  recent  progress  in  our  analysis  of  the  EAW  communication 
scheme.  Sec.  II,  and  the  ReComm  scheme.  Sec.  III. 

II.  ELECTRON  ACOUSTIC  WAVE  COMMUNICATION 

One  scheme  for  commnnication  through  a  plasma  sheath  is  based  on  the  possibility  of 
propagating  EAWs  with  freqnencies  well  below  the  local  plasma  frequency  (~  9  GHz)  in  a 
two-electron-temperature  plasma.  A  small  hot  electron  population  is  produced  by  injection 
of  an  energetic  electron  beam.  Over  a  certain  range  of  wavelengths,  the  hot-electron  com¬ 
ponent  of  the  plasma  can  short  out  the  electric  helds  that  produce  plasma  oscillations  in 
the  cold  plasma,  thereby  redncing  the  freqnency  of  these  oscillations  below  the  cold  plasma 
frequency.  It  may  be  possible  to  excite  these  electrostatic  oscillations  with  an  antenna  and 
have  them  conple  to  electromagnetic  RE  oscillations  ontside  the  sheath.  Related  research 
was  carried  out  in  the  CHARGE  2B  ionospheric  rocket  experiments  [18-20].  For  sustained 
hypersonic  flight,  a  comprehensive  analysis  is  reqnired  to  estimate  the  feasibility  of  this 
concept. 

Investigations  into  the  EAW  communication  scheme  have  yielded  a  dispersion  analysis 
indicating  that  EAW  modes  can  be  generated  and  propagate  in  the  plasma  layer.  Further¬ 
more,  we  have  demonstrated  that  these  modes  can  be  conpled  to  electromagnetic  waves  at 
the  plasma  layer  boundary  for  the  assumption  of  a  sharp  density  bonndary.  We  presented 
a  detailed  analytic  analysis  of  the  wave  transformation  coefficients  at  this  bonndary  layer 
in  Ref.  [11].  We  have  begun  developing  a  nnmerical  simulation  model  of  EAW  wave  gen¬ 
eration  and  propagation  for  testing  the  characteristics  of  these  waves.  Below  we  describe  a 
particle-in-cell  (PIC)  simulation  model  presently  under  development  for  analyzing  the  EAW 
commnnication  scheme. 
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A.  Simulation  model 


The  initial  PIC  simulations  described  here  are  highly  idealized  ID  representations  of  an 
electron  beam  propagating  through  a  two-temperature  fully  ionized,  collisionless  plasma. 
The  electron  beam  is  the  free  energy  source  that  drives  unstable  modes  in  the  plasma, 
including  EAWs.  In  the  near  future,  we  will  extend  this  model  to  2D  to  include  the  effects 
of  electron-neutral  collisions  and  hnite  plasma-sheath  thickness  to  examine  the  coupling  of 
EM  waves  to  EAWs  for  a  variety  of  plasma  prohles.  We  anticipate  that  these  2D  simulations 
will  be  large  scale,  so  the  ID  simulations  presented  here  are  a  logical  hrst  step  in  the  overall 
analysis  of  the  EAW  communication  scheme  to  test  numerical  constraints. 

The  ID  simulation  model  assumes  a  periodic  system  of  length  L,  with  a  preformed  plasma 
and  inter-penetrated  electron  beam.  The  electron  beam  has  a  small  but  hnite  thermal  spread 
Tfe  which  is  less  than  the  beam  speed  Uf,.  The  electron  beam  number  density,  rib,  is  smaller 
than  the  background  electron  number  density,  rip.  The  background  plasma  is  composed  of 
hot  and  cold  electron  populations  with  temperatures,  Th  and  Tc,  and  number  densities, 

np  =  nh  +  nc  =  ni-  rib,  (1) 

where  n*  is  the  background  ion  density.  For  the  initial  simulations  described  here,  we 
assume  a  cold,  neutralizing,  immobile  ion  background  population  with  charge  state  Z  =  1. 
Extension  to  a  thermal  ion  background  population  does  not  affect  the  results  of  these  short 
time  duration  simulations  but  can  be  trivially  included  in  future  studies. 

The  LsP  simulation  code  is  used  in  ID  with  the  explicit  EM  held  solver.  We  note  that 
these  simulations  can  also  be  carried  out  in  an  electrostatic  (ES)  mode,  but  our  goal  is  the 
development  of  the  general  2D  model  that  includes  both  ES  and  EM  wave  propagation,  and 
therefore  we  carry  out  the  ID  modeling  using  the  EM  held  solver.  These  ID  simulations  are 
carried  out  in  parallel  on  4  to  16  processors  and  extension  to  2D  will  likely  require  additional 
processors. 

The  simulations  are  initialized  with  400  particles  per  species  per  grid  cell  to  provide 
good  velocity  space  resolution  of  evolving  particle  energy  distribution  functions  throughout. 
We  benchmarked  our  simulations  against  recent  work  by  Q.  Lu  et  al.  [21]  where  similar 
ID  ES  PIC  simulations  were  carried  out  to  test  electron  beam  generated  EAW  dispersion 
characteristics.  Good  agreement  between  our  EM  simulations  and  these  previously  published 
ES  simulations  was  found  for  all  cases. 
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B.  Sample  simulation  results 


Our  initial  tests  have  focused  on  one  of  the  parameter  sets  from  our  last  report  (e.g., 
case  1).  We  consider  a  two-temperature  plasma  with  Tc  =  1  eV,  Th  =  20  eV,  and  number 
densities  Uc  =  8.6  x  10^^  cm“^  and  Uh  =  5  x  10^°  cm“^.  The  electron  beam  parameters 
are  Tf,  =  0.1  eV,  rib  =  9  x  10^°  cm“^,  and  Vb  =  0.0267c,  where  c  is  the  speed  of  light.  The 
background  ion  density  n*  is  then  1  x  10^^  cm“^.  The  simulation  length  L  is  1.862  cm, 
divided  into  2048  equal  size  grid  cells,  giving  a  grid  size  of  Ax  ~  9.09  x  10“^  cm.  The 
constraint  on  numerical  grid  heating  [22] 


Ax  <  TrXoe, 


(2) 


where  is  defined  as  A^e  —  743  (Te/ue)^'^^,  is  well  satished  for  Te  >  1  and  Ue  =  1  x  10^^  cm 
The  time  step  for  these  explicit  simulations  is  At  ~  2.97  x  10“^^  s,  which  resolves  the  electron 
plasma  frequency 

UpeAt  <  ^,  (3) 


where  Upe  ~  5.64  x  10"^  is  the  electron  plasma  frequency.  This  time  step  also  resolves 

the  speed  of  light  in  vacuum  across  the  individual  cells. 

The  initial  particle  velocity  distributions  can  be  seen  in  Fig.  la),  which  is  a  plot  of  the 
x-Vx  phase  space  at  t  =  0.  Only  the  first  1-cm  of  the  simulation  is  shown  to  more  clearly 
illustrate  the  beam-plasma  interaction.  The  beam  electrons  are  shown  in  green,  the  hot 
electrons  are  shown  in  red,  and  the  cold  electrons  are  shown  in  blue.  The  cold  electrons 
are  plotted  first  in  this  figure,  so  over-plotting  of  the  hot  electrons  largely  obscures  the  cold 
electrons.  The  small  thermal  spread  of  the  beam  electrons  is  clear  when  compared  with  the 
hot  electrons  at  t  =  0.  The  temporal  evolution  of  the  electron  phase  space  grows  initially 
according  the  beam-plasma  instability,  which  has  a  maximum  unstable  wavenumber  given 
approximately  by 


^max  —  ? 

Vb 


(4) 


which  for  our  parameters  is  roughly  kmax  —  70  rad/cm,  or  0.09  cm.  At  time  Upet  ~  23 
in  frame  (c),  the  electron  beam  has  developed  a  sinusoidal  modulation  with  wavelength  of 
about  0.08  cm,  which  is  quite  close  the  approximate  calculation  of  the  maximum  unstable 
wavenumber.  In  frames  (d)  through  (f)  of  Fig.  1,  the  nonlinear  phase  of  the  beam-plasma 
instability  continues  to  grow,  illustrating  electron  trapping  which  leads  to  the  saturation  of 
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FIG.  1:  Phase  space  plots  of  beam  (green),  hot  plasma  (red)  and  cold  plasma  (blue)  electrons  at 
different  times  from  a  ID  PIC  simulation  of  EAW  generation. 

the  instability.  After  approximately  ojpet  ~  100,  the  beam  electrons  have  largely  thermalized 
with  the  hot  and  cold  electron  popnlations. 

The  late  time  stage  of  the  beam-plasma  interaction  as  modeled  here  is  not  exactly  repre¬ 
sentative  of  the  hypersonic  commnnication  scheme,  since  the  electron  beam  injected  into  the 
plasma  sheath  will  essentially  pass  throngh  the  3-5  cm  thick  plasma  sheath  and  propagate 
away  from  the  vehicle.  Nevertheless,  the  ID  simnlation  provides  a  model  for  self-consistently 
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kilovolts/cm 


X  (cm) 

FIG.  2:  Logio  contours  of  \Ex\  as  a  function  of  time  (ujpet)  and  space  (x)  from  the  ID  PIC 
simulation. 

generating  EAWs  under  a  variety  of  ambient  plasma  conditions. 

The  generation  of  the  EAWs  can  be  examined  by  plotting  the  electric  held  from  the 
simulations  as  a  function  of  space  and  time.  Figure  2  plots  the  magnitude  of  the  electric 
held  in  the  simulation  as  a  function  of  time  {upet)  and  space  (x).  From  this  data  we  examine 
frequencies  and  wavelengths  for  EAW  modes. 

Figure  3  plots  the  Fourier  transform  of  the  total  magnetic  held  energy  as  a  function  of 
time  in  the  simulation  volume.  The  two  dominant  features  of  the  frequency  spectrum  are  the 
electron  plasma  frequency  {cOpe  ~  9  x  10^  MHz)  and  the  EAWs.  From  the  EAW  dispersion 
analysis  presented  in  [11],  the  real  frequencies  of  the  EAWs  are  roughly  centered  around 
^EAw  ~  0.2ujpc  and  over  wavenumbers  ranging  from  kXnh  —  0.1  to  0.3  for  these  parameters. 
The  simulation  presented  here  was  only  run  out  to  about  250  plasma  periods,  which  gives 
only  about  8  EAW  periods. 


frequency  (MHz) 

RMS  Amplitude  =  1.17430e-007  in  2  bins 


FIG.  3:  Fourier  transform  of  the  total  electric  field  energy  in  the  simulation  volume  as  a  function 
of  time.  The  electron  plasma  frequency  and  the  EAW  regions  are  indicated. 

III.  ANALYSIS  OF  RECOMM  SCHEME:  CROSSED-FIELD  DIODE 

The  ReComm  scheme  utilizes  externally  generated  crossed  electric  and  magnetic  helds  to 
induce  charged-particle  drifts.  This  spatially  redistributes  the  plasma  and  enables  propaga¬ 
tion  of  electromagnetic  waves  through  the  plasma  layer.  An  electric  held  is  applied  across 
a  diode  placed  in  the  plasma  sheath.  In  previous  ID  simulations  [14]  we  demonstrated  that 
the  time  required  to  signihcantly  alter  the  plasma  density  prohle  in  the  diode  corresponds 
roughly  to  the  ambipolar  diffusion  time.  When  a  voltage  is  applied,  the  plasma  center-of- 
mass  shifts  from  the  center  of  the  diode.  But  whether  that  shift  is  toward  the  cathode  or 
anode  depends  on  the  relative  mobility  of  the  plasma  ions  and  electrons  in  the  non-neutral 
sheaths  at  the  electrodes. 

A  simplihed  picture  of  the  ReComm  scheme  can  be  seen  in  Fig.  4.  The  antenna  aperture 
is  placed  in  between  a  pair  of  biased  electrodes.  An  electromagnet  provides  a  magnetic  field 
normal  to  the  antenna  aperture  throughout  the  region  between  the  electrodes.  We  have 
considered  the  “plasma-optic”  regime  described  by  Keidar  et  al.  [23]  in  which  electrons  are 
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magnetized  and  ions  are  unmagnetized.  The  separation  distance  between  the  plates,  L,  is 
assumed  to  be  large  compared  to  the  electron  Larmor  radius,  but  small  with  respect  to  the 
ion  Larmor  radius.  In  this  case,  if  E  is  the  electrostatic  held  between  the  electrodes,  the 
electrons  are  line-tied  and  can  only  E  x  B  drift  parallel  to  the  electrodes.  In  principle  [23,  24] , 
this  prevents  the  electrons  from  screening  out  the  electric  held  in  the  bulk  of  the  plasma 
between  the  electrodes.  The  ions,  being  unmagnetized,  are  then  accelerated  ballistically 
in  the  unscreened  electric  held.  This  simple  description  assumes  that  electrons  and  ions 
are  collisionless.  If  there  is  a  hnite  collision  frequency,  the  electrons  can  dihuse  across  held 
lines.  Similarly  ion  acceleration  across  the  electrodes  will  be  limited  by  ion  collisionality. 
Including  both  ion  and  electron  collisionality  in  a  ID  huid  theory  analysis,  Keidar  et  al.  [23] 
argue  that  a  steady-state  spatial  distribution  results  in  which  the  plasma  density  is  reduced 
near  the  cathode  side.  The  amount  of  density  reduction  will  be  a  function  of  the  plasma 
properties,  electrode  voltage,  and  magnetic  held  strength.  The  aperture  of  the  antenna  used 
for  communications  would  then  be  placed  directly  below  the  reduced  density  region,  where 
the  plasma  frequency  will  be  lower.  A  communication  frequency  of  1  GHz  corresponds  to  a 
free-space  wavelength  of  30  cm.  The  antenna  aperture  will  be  on  the  order  of  the  wavelength, 
which  means,  as  a  rough  estimate,  we  can  assume  that  the  electrode  separation  L  must  be 
on  the  order  of  a  meter. 

In  Keidar’s  huid  theory  analysis  he  neglects  the  time- variation  in  the  continuity  equations, 
that  is  he  considers  timescales  short  compared  to  the  ambipolar  dihusion  time,  and  hnds  a 
steady-state  solution  by  retaining  convection  in  the  ion  momentum  equation.  In  our  PIC 
simulations  we  have  observed  no  signs  of  such  a  steady-state  density  reduction  on  short- 
timescales.  Rather  we  found  that  on  very  short  (ion)  timescales  {t  ~  Upi)  non- neutral  Debye 
sheaths  are  established  on  the  electrode  boundaries.  This  is  followed  by  dihusive  ehects 
which  occur  on  a  relatively  long  ambi-polar  dihusion  timescale  td.  Moreover  we  were  able 
to  qualitatively  explain  the  simulation  density  and  potential  prohles  with  a  two-huid  model 
in  the  dihusion  approximation,  neglecting  ion  convection. 

The  LsP  PIC  simulations  consist  of  a  ID  grid  with  conductors  at  the  boundaries.  A 
constant  voltage  diherence  (after  a  short  temporal  ramp)  of  V  is  applied  between  the  anode 
(x  =  0)  and  cathode  {x  =  L).  Table  I  shows  the  physical  parameters  used  in  what  we  will 
refer  to  as  the  “baseline”  ReComm  simulation.  The  electron-collision  frequency  value  used 
in  the  baseline  simulation  corresponds  to  a  value  of  Uen  —  0.02upe,  which  is  in  accord  with 
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FIG.  4:  ReComm  communication  through  plasma  layer  using  applied  electric  and  magnetic  fields. 
Adapted  from  Ref.  [25]. 

estimates  of  the  value  in  the  hypersonic  plasma  sheath  [6,  13].  Since  we  have  not  found  any 
experimental  data  or  estimates  for  the  ion-neutral  collision  frequency,  we  have  chosen  to  set 
^in  =  where  Upi  is  the  ion  plasma  frequency.  For  an  applied  magnetic  held  of  2000  G, 
Venl'^ce  ~  0.01.  By  coutrast  the  ions  are  in  the  opposite  regime:  VinloJd  ~  100.  Since  the 
electron  Larmor  radius  is  r^e  —  0.04  cm,  the  electrons  will  be  line-tied  at  the  nominal  plate 
separation  of  L  =  1  cm. 

By  choosing  a  plate  separation  of  1  cm,  we  can  choose  a  cell  size  of  Ax  =  0.04  cm, 
which  is  on  the  order  of  the  electron  cyclotron  radius,  without  requiring  a  prohibitively 
large  number  of  computational  grid  cells.  Several  hundred  plasma  particles  per  cell  were 
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Plate  separation 

L  (cm) 

1.0 

Applied  voltage 

F(V) 

100 

Applied  magnetic  field 

B{G) 

2000 

Initial  plasma  density 

n  (cm“^) 

10^1 

Initial  Electron  temperature 

Te  (eV) 

10.0 

Initial  Ion  temperature 

Ti  (eV) 

0.1 

Ion-neutral  collision  frequency 

Uin  (ns“^) 

0.066 

Electron-neutral  collision  frequency 

Uen  (ns“^) 

0.329 

TABLE  I:  Physical  Parameters  for  the  baseline  ID  ReComm  simulation.  A  spatially  uniform  singly 
ionized  Ar  plasma  is  placed  between  plate  electrodes. 

used  in  these  simulations  to  provide  a  good  statistical  representation  of  the  electron  and  ion 
energy  distribution  functions.  A  time  step  of  At  =  0.025  ns  gives  UceAt  ~  2ujpeAt  ~  0.9, 
ensuring  numerical  stability.  Simulations  with  a  signihcantly  larger  time  step  would  require 
the  use  of  implicit  algorithms  (Ref.  [26]). 

Figure  5  shows  the  results  of  a  series  of  ID  simulations  in  which  the  electrode  separation 
L  is  varied.  All  other  parameters  are  the  same  as  those  given  in  Table  I.  All  the  simulations 
neglect  the  flow  velocity  of  the  plasma.  In  a  hypersonic  vehicle  the  flow  velocity  could  be 
as  large  as  Mach  10,  with  a  sound  speed  on  the  order  of  300  m/s.  For  the  ReComm  scheme 
to  be  useful  in  practical  applications,  the  diffusion  time  td  must  be  short  compared  to  the 
plasma  £11  time  rgjj  =  L/Cg,  where  Cg  is  the  local  sound  speed.  The  fill  time  is  just  the  time 
required  for  a  plasma  moving  at  the  flow  velocity  to  fill  the  region  between  the  electrodes. 
If  Tgjj  is  not  large  compared  to  td  plasma  will  flow  into  the  electrode  region  faster  than  it 
can  be  diffused  away. 

In  Fig.  5  we  plot  the  diffusion  time  as  a  function  as  a  function  of  L.  There  are  two 
separate  diffusion  times  plotted  in  the  figure,  which  are  calculated  by  the  formulas: 

Q{t)  =  Const  X  exp  ( — - 

V^fit 

Qi^e)  =  Q{0)e~\  (5) 

where  Q{t)  is  the  total  ion  charge  in  the  simulation  (neglecting  plasma  refill).  The  figure 
shows  that  for  small  L,  tb  scales  roughly  as  while  at  larger  values  of  L  the  scaling  is 
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L  (cm) 


FIG.  5:  Results  from  a  series  of  ID  ReComm  simulations.  The  simulation  parameters  are  the  same 
as  those  given  in  Table  I,  except  that  the  plate  separation  L  is  varied.  Diffusion  timescales  Tf  and 
Te  are  plotted  as  a  function  of  L. 

closer  to  linear.  Both  the  linear  and  quadratic  fits  are  extended  out  to  larger  values  of  L  in 
the  figure  and  are  shown  as  dashed  lines.  The  vertical  black  line  indicates  L  =  1  ni,  which 
is  of  the  order  of  the  desired  plate  separation,  while  the  horizontal  line  shows  the  plasma  £11 
time  for  a  Mach  10  plasma  for  the  same  separation.  From  the  figure  we  see  that  even  for 
the  more  favorable  linear  scaling  the  diffusion  time  at  L  =  1  m  exceeds  the  plasma  fill  time 
by  an  order  of  magnitude. 

We  have  also  performed  a  series  of  simulations  with  varying  values  of  applied  voltage  (see 
Figure  8  of  [14]).  In  this  series  of  simulations  we  found  (with  a  fixed  value  L  =  1  cm)  that 
both  and  Tg  scale  roughly  as  V~^.  If  we  assume  that  this  scaling  with  voltage  holds  up 
at  larger  values  of  L,  we  could  obtain  td  <C  igjj  by  increasing  the  voltage  from  the  nominal 
value  of  100  V  by  two  to  three  orders  of  magnitude. 

In  the  future  we  plan  to  do  simulations  to  demonstrate  that  the  ReComm  scheme  can  be 
successfully  employed  at  L  ~  1  m  with  a  reasonable  applied  voltage  .  As  mentioned  above, 
without  converting  to  implicit  algorithms  these  will  be  considerably  longer  simulations  than 
those  previously  performed  with  short  values  of  L.  We  also  re-emphasize  our  conclusion 
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from  Ref.  [14]  that  the  shape  of  the  density  prohle  depends  on  the  values  of  electron  and 
ion  collision  frequencies  compared  to  cyclotron  frequencies.  Ultimately,  we  will  need  more 
accurate  values  of  Ven  and  Ven  in  the  sheath  to  accurately  predict  density  prohles. 

If  it  is  established  that  the  ReComm  scheme  can  be  usefully  employed  in  ID,  we  would 
like  to  go  on  to  consider  multidimensional  effects.  In  the  ID  simulations  we  have  assumed 
that  the  electric  held  throughout  the  plasma  is  only  a  function  of  x  (see  Fig.  4).  This  neglects 
the  effect  of  the  hnite  electrode  and  plasma  thicknesses  as  well  as  fringe  helds  and  boundary 
effects  at  the  vehicle  surface.  For  example,  if  the  electrodes  are  placed  below  the  surface  of 
the  vehicle  behind  a  dielectric  radome,  will  the  electric  held  penetrate  a  sufficient  depth  into 
the  plasma  layer?  Such  multi-dimensional  ehects  can  be  investigated  by  massively  parallel 
2-D  LsP  simulations.  We  have  also  previously  demonstrated  [12]  that  the  presence  of  a  large 
magnetic  held  alone  can  considerably  increase  the  transparency  of  the  plasma  sheath  even 
in  the  absence  of  an  applied  voltage.  It  would  be  interesting  to  consider  the  combined  ehect 
of  the  “magnetic  window”  and  plasma  redistribution  due  to  crossed-helds. 

IV.  DISCUSSION  AND  SUMMARY 

The  analysis  of  the  conversion  between  EM  and  EAWs  has  been  carried  out  for  a  sharp 
plasma  boundary,  and  the  results  demonstrate  reasonable  power  for  the  transmitted  EM 
waves.  We  are  now  developing  an  EM  PIC  simulation  model  to  assess  the  generation  and 
propagation  of  the  EAWs  and  their  transformation  to  EM  waves  at  the  plasma-vacuum 
interface. 

Our  initial  analysis  of  the  ReComm  scheme  [23,  27,  28],  where  applied  electric  and  mag¬ 
netic  helds  are  used  to  open  a  frequency-space  window  for  EM  wave  propagation  through 
the  plasma  sheath,  is  presented  in  Ref.  [14].  We  established  that  the  ion  distribution  in 
the  cross-held  dipole  (in  the  plasma-optic  regime)  is  altered  signihcantly  only  on  dihusion 
timescales,  for  which  r  oc  jV ^  where  7  ~  1  —  2.  Eor  plasma  parameters  believed  to  be 
consistent  with  the  radio  blackout  regime,  r  is  of  the  order  of  tens  of  ps,  when  L  on  the 
order  of  a  few  cm,  and  V  is  on  the  order  of  100  V.  For  a  realistic  ReComm  device  allowing 
communication  in  the  L-band  of  the  radio  spectrum,  the  diode  length  L  would  have  to  be 
scaled  up  to  something  on  the  order  of  1  m.  The  magnitude  of  the  applied  voltage  V  is, 
of  course,  limited  by  practical  considerations.  The  utility  of  the  ReComm  scheme  depends 
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upon  being  able  to  alter  the  plasma  density  profile  on  a  timescale  which  is  short  compared 
to  the  time  it  takes  the  hypersonically  flowing  plasma  to  pass  by  the  diode. 

We  have  shown  that  an  ambipolar  diffusion  theory  of  the  cross-held  diode  can  describe 
qualitatively  the  results  of  the  simulation.  However,  plasma  energy  transport  and  charge 
separation  effects,  neglected  in  the  theory,  play  an  important  role  in  determining  the  details 
of  the  plasma  evolution.  We  have  shown  in  our  baseline  simulation  that  ion  temperatures 
can  reach  values  of  several  tens  of  eV  on  long  timescales  t  ~  r.  This  snggests  that  Ohmic 
heating  of  the  plasma  may  lead  to  further  ionization.  To  examine  snch  effects  would  require 
replacing  the  simple  Drnde  collision  model  used  in  this  report  with  a  more  realistic  chemistry 
model.  Although  it  is  probably  more  important  to  establish  initially  that  r  can  be  made 
small  enough  for  practical  use. 

Another  concern  is  that  the  direction  in  which  the  plasma  is  shifted  by  the  applied  held 
depends  sensitively  on  the  relative  mobility  of  the  electrons  and  ions.  To  determine  the 
relative  mobility,  it  is  necessary  to  have  reliable  information  about  the  electron  and  ion 
temperatnres  and  collision  freqnencies  in  the  sheath  surronnding  the  vehicle.  Onr  survey  of 
the  literature  suggests  that  the  existence  of  this  information  is  still  incomplete. 

We  note  that  the  ReComm  scheme  is  based  only  on  altering  the  local  plasma  density 
in  the  vicinity  of  the  antenna.  To  do  this  reqnires  both  an  applied  magnetic  held  and  a 
high  voltage  diode.  Moreover  the  ehect  may  occnr  only  on  prohibitively  long  timescales.  By 
contrast  the  magnetic  window  approach  [29],  which  ntilizes  the  magnet  bnt  does  not  reqnire 
the  diode,  does  not  require  any  slow  bulk  motion  of  the  plasma,  but  instead  alters  the 
transparency  of  the  existing  plasma  to  the  RF  signal.  But  the  ehectiveness  of  this  alternate 
scheme  depends  strongly  on  the  plasma  collisionality  [13],  which,  as  noted  above,  is  not  well 
known. 

In  all  communication  scenarios  under  consideration,  the  computational  models  will  be 
extended  to  include,  with  increasing  levels  of  complexity,  realistic  gas  and  plasma  properties 
associated  with  hypersonic  hight  [30-32].  The  goal  of  this  integrated  theoretical  analysis  is 
the  formulation  of  “operational-windows”  in  which  physical  constraints  associated  with  each 
commnnication  scheme  can  be  parameterized.  This  analysis  should  signihcantly  impact  the 
down-select  of  these  schemes  for  eventual  deployment  on  hypersonic  vehicles. 
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